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1 Introduction 

It is usual to consider Hamiltonian systems in E 3 of the form H = + /(r) + \L, where L = 
2 (Pg + iE7gp' P^)"l~^ r (^i VO) an d i r ,@,4>) are spherical coordinates. The examples are innumerable 
and we just cite the n-body systems of points in one dimension pQ whose Hamiltonian can 
be written in a similar way for some L on § n . In this paper we restrict our attention to 
Hamiltonians of the form H = \p\ + f(u) + a(u)L in n + 1 dimensions that are extensions of 
a n-dimensional Hamiltonian L(q l ,pi), by calling u the (n + l)-th coordinate and p u its canonical 
conjugate momentum. We will call them extended Hamiltonians for the given L. We show how 
to construct a non trivial first integral of H in the form U m (G), m £ N*, where U is a differential 
operator built from the Hamiltonian flow of L, and G is a suitable function. If L is a natural 
Hamiltonian and G is independent of the momenta, the first integral U m {G) is a m-th degree 
polynomial in the momenta (p u ,Pi)- If the Hamiltonian L is integrable, then automatically H 
has n + 2 functionally independent first integrals and is therefore minimally superintegrable. 
By requiring certain natural integrability conditions, we show that a necessary condition for the 
construction is the constant curvature equal to mc of the metric tensor of L, where c is a real 
constant. In few words, for any non null positive integer m, given a n-dimensional Riemannian 
or pseudo- Riemannian manifold of constant curvature equal to mc, underlying the geodesic part 

*This paper is a contribution to the Special Issue "Symmetry, Separation, Super-integrability and Special 
Functions (S 4 )". The full collection is available at |http://www.emis.de/journals/SIGMA/S4.html| 
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of the Hamiltonian L, our procedure can construct a natural Hamiltonian H inn + 1 dimensions 
with an additional first integral U m (G) which is polynomial of degree m in the momenta. If L 
includes a non constant scalar potential V, compatibility conditions on V are given in order that 
the procedure may work, otherwise, all the compatible V can be determined from a first-order 
PDE. This approach is motivated by the results obtained in [5], where L was a one-dimensional 
Hamiltonian. In that paper, we obtained for L = \p^ + V(ip) and a = u~ 2 the additional first 
integral 

in 

cos{mip + ^o) 

with V(ip) = sin ^( m *0 + ^, o ) i an d G(i/}) = cos(rmp + ipo). The resulting Hamiltonian H can be 
interpreted as generalization of the three-particle Calogero and Wolfes systems and is maximally 
superintegrable jl]. It represents also an instance of the more general Tramblay-Turbiner- 
Winternitz systems \15\ [T6] whose superintegrability has been proved by another way in . 
In recent years much research has been done about superintegrable natural Hamiltonians with 
polynomial first integrals of high degree in the momenta, both classical an quantum. Many 
results concern systems which admit separable Hamilton-Jacobi equations and the separability 
of those systems is explicitly employed to build the high-degree first integrals, see for example [HJ 
El [HI ESI [17] . The procedure presented here to build the additional first integral does not assume 
separability anywhere. However, the requirement of working with constant curvature manifolds 
allows to consider the extensions of separable Hamiltonian systems, because constant curvature 
manifolds are the most natural environment for those systems (see for example [7]). It is not by 
accident that some of the examples of superintegrable Hamiltonians that we provide here in the 
section of two and three-dimensional systems reveal themselves to be extensions of separable 
systems. 



U m (G) 



Pu + 



mu 



d dV d 

^ dtp dip dp^, 



2 Main results 

Let us consider a Poisson manifold M and a one-dimensional manifold N. For any Hamiltonian 
function L G J~(M) with Hamiltonian vector field Xl, we consider its extension on M = 
T*N x M given by the Hamiltonian 

H=±p 2 u + a(u)L + /(«), (1) 
where (p u ,u) are canonical coordinates on T*N. The Hamiltonian flow of ([T]) is 
X H = Pu— - (a'L + f)g^- + uX L , 

where prime denotes total derivative w.r.t. the corresponding variable. 

It is immediate to see that any first integral of L is also a constant of motion of H, when 
considered function on M. 

We want to determine on L, a and / necessary and sufficient conditions for the existence of 
two functions 7 G F(N) and G € F{M) such that, given the differential operator 

U = Pu + 1 {u)X L , (2) 

the function F obtained applying m times U to G 

F = U m {G) 



is an additional first integral for H. 
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Remark 1. We observe that the operator ([2]) is tensorial with respect to transformations acting 
on M only i.e. preserving the form of the extended Hamiltonian (TT]). Moreover U is an injective 
linear operator on the space of functions polynomial in the momentum p u : indeed if F £ F(M) 
is polynomial in p u then U(F) = if and only if F = because U(F) is a polynomial of degree 
higher than F. 

Proposition 1. For any L, a, f , we have that X[{U ni {G) = for a function G(q l ,pi) if and 
only if 

{mi + a)X L {G) = 0, (3) 
a-fXl (G) - m(a'L + f)G = 0. (4) 

Proof. Being 

[X H ,U]= Pu j'X L -La' -f, 

we have £/],{/] = 0. Hence, we can apply the following formula for the power of a differen- 

tial operator: 

X H U rn = mU m -\X H , U] + U m X H = U m -\m[X H , U] + UX H ), 
Since we have 



2 d 

UX H =Puq^ +Pu 



OPu OU 



+ a 1 X 2 L - 1 { a 'L + f)X~, 

OPu 



and, in particular, for any G(q l ,pi), 

UX H {G) = a(p u X L (G) + 7 X 2 L (G)), 
thus, we get for XHU m (G) the following expression 

U m - 1 (p u {m-/ + a)X L (G) + ajX 2 L (G) - m(a'L + f')G). 

Since the operator U is injective on the space of functions polynomial in p u we have the thesis. ■ 

The conditions ([3]) and (H} are quite difficult to handle in the general case, their analysis can 
be further pursued by adding some extra hypotheses on the manifold M and on the Hamilto- 
nian L. 

Theorem 1. Let Q be a n- dimensional (pseudo-) Riemannian manifold with metric tensor g. 
The natural Hamiltonian L = hg^PiPj + V(q l ) on M = T*Q with canonical coordinates 
admits an extension H in the form JT]) with a first integral F = U m (G) with U given by ([2]) 
and G{q l ), if and only if the following conditions hold: 

1) the functions G and V satisfy 

H(G) + mcgG = 0, ceR, (5) 
W • VG - 2m{cV + L )G = 0, Lq G R, (6) 

where H(G)y = VjVjG is the Hessian tensor of G. 
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2) for c = the extended Hamiltonian H is 



H = 




for c the extended Hamiltonian H is 



h = \pI + 



m{cL + Lq) 



+ W 



Sl{cu + u ) 



with k, u , Vq, W ,A e R, A ^= 0, B = mL A 2 and 



' sin y/Tix 



K > 



S K {x) 



= < 



X, 

sinh 



K = 




K < 0. 



For the equations © it is natural to write integrability conditions. It follows (Lemma [3|) that 
the maximal dimension of the space of solutions of equation ([5]) alone is n + 1 and is achieved 
only in the constant curvature case. We call complete the solutions G of ([5]) satisfying these 
integrability conditions. Therefore, we can restate Theorem [1] as follows. 

Theorem 2. Let Q be a n-dimensional (pseudo-) Riemannian manifold with metric tensor g. 
The natural Hamiltonian L = hg^PiPj + V(q l ) on M = T*Q with canonical coordinates {pi,q l ) 
admits an extension H in the form ([T]) with a first integral F = U m (G) with U given by ([2]) 
and G(q l ) a complete solution of ©, if and only if Q is a (pseudo) -Riemannian manifold with 
constant curvature K = mc and ([6]) holds. 

Clearly, in the case of Theorem [2] the extended Hamiltonian H is again in the form given by 
item 2 of Theorem [TJ 

Depending on the form of V, the space of common solutions of the two equations ([5]) and © 
can be of lower dimension but, in any case, the maximal number of independent solutions is 
obtained in the constant curvature case. Since ([5]) and © are linear differential conditions on G, 
their common solutions form a linear space of dimension lower or equal to n + 1 (dimension of 
the space solution of © alone) parametrized by up to n + 1 real parameters (oi,. . . ,a n+ i); 
the linear injective operator U k , k G N*, maps this space in a new linear space, with the same 
dimension, spanned by the functions U k (G). These functions are first integrals of H if and 
only if k = m. It is possible that functions U m (G) corresponding to different choices of the 
parameters (oj) are simultaneously first integrals of H and functionally independent even for 
a non-trivial potential V . The analysis of this interesting case and of the relations among 
first integrals in general will be done in some future paper. The condition ([6]) can be read as 
a condition on the potentials V, whose solutions form a functional space of potentials depending 
on (at most) n + 1 parameters and admitting at least one non zero G such that F = U m (G) is 
a first integral. Examples of potentials V depending on some of the parameters (oj) are given in 
the next section. The presence of parameters into the potentials V of superintegrable systems 
can allow the application of the Stackel transform or, more generally, of the coupling constant 
metamorphosis to obtain new superintegrable systems, see for example |10j and |14| . 

Remark 2. Equation ([5]) could admit in particular cases non- vanishing solutions also on man- 
ifold with non-constant curvature, provided they depend on less than n + 1 parameters and, as 
appears from some example, do not depend on some of the variables (q l ). A detailed analysis 
of this case is in progress. 
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The proofs of the theorems follow from the next three lemmas. 

Lemma 1. If L is a natural Hamiltonian and G a function of (q l ) we have XH~U m {G) = - 
with H and U defined as in ([T]) and ^ - if and only if the functions a, j, f , G and V satisfy 
the following differential conditions 

a = —rwy', (7) 

(y + cr > y = o, (8) 

V i V j G = -mcGg ij , (9) 
ViVVG = 2m(cV + L )G, (10) 
/ = mLol 2 + f , (11) 

where c, Lq and fo are arbitrary constants. 

Proof. By using the (pseudo-)Riemannian structure of Q and covariant derivatives we have 
X L =p i V i -V i V^-, 

OPi 

X\ = pffjV'V' - V l VV i - 2 Pj ViVV j ^- - ftV'V//- + ViVVjV- 



dpi dpj dpidpj 

Hence, if G does not depend on the momenta, we have 

X L (G) = Pl V l G, X 2 L (G) = p^V^G - ViVV'G. 

By substituting in ([3]) and excluding the trivial case of constant G we obtain (J7]) and the two 
conditions 

777 

orfVtVjG = -a' 9ij G, (12) 

ajX7 i VV i G + m(a'V + f')G = 0. (13) 

From (|12p . by separating terms in u from the other ones we get (|9|) and 

a' = -2ca7, (14) 

where c is a real constant. By substituting ([7]) in the relation (|14p we obtain 7" + 2077' that 
is dU). By p^|) the equation (fT3|) is equivalent to 



aj(VV ■ VG - 2mcVG) + mf'G = 0, 
which by separating terms in u from the others and substituting ((7|) splits into (|10p and 

/' = -2L «7 = 2mL 77', 
with Lq constant, that integrated gives (jlip . ■ 

Remark 3. The addition of a term 4>(u) to U does not induce weaker constraints on the func- 
tion f(u) in dl]): a straightforward calculation shows that (j>(u) must necessarily vanish. 

Condition ((HJ) of Lemma [T] determines the possible forms of the extension. The function 7 
and the related extended Hamiltonian are given by the next lemma and are expressed by using 
the "tagged" trigonometric functions 

' sin ^/kx 

cos ^JUx^ k > 0, 
1, K = 0, S K (x) = 

cosh a/|k|x, k < 0, 



K > 0, 



IK 

X, K = 0, 

sinh 

-, K < 0, 
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T K (x) 



S K (x) 
C K (x)' 



CT K (x) 



C K (x) 
S K (x) 



already employed in [13|. O [12]. These function share almost all the properties of the stan- 
dard trigonometric and hyperbolic functions, conveniently modified with the tag k, as some 
straightforward calculations will show. 

Lemma 2. The equation 7" + 2077' = for c = reduces to 7" = 0, its solution is 7 = 
— A(u + «o) with uq and A / constants. The related extended Hamiltonian is 

H = \pl + mA(L + Vq) + B(u + u ) 2 , 



m 



LqA 2 , i.e. the sum of L and the Hamiltonian of a harmonic 



with uq, Vb, A e R, A ^ 0, B 
oscillator. 

The equation 7" + 2077' = for c 7^ is equivalent to 7' = — c(7 2 + k) with k a constant, its 
solution is: 



1 



1 



T K (cu + u ) 
and the related extended Hamiltonian is: 



where Wq € 



1 2 m(cL + Lq) 



52 (cm + n ) 



Proof. The first part of the proposition is trivial: equations ([7|) and (jlip give a = mA and 
/ = mLQA 2 (u + uq) 2 + /o- The constant A must be assumed different from zero to avoid the 
vanishing of a. The given form of H is obtained by setting B = uiLqA 2 and Vo = tjjjr. 

The second part immediately follows from the fact that if CT~ 1 (x) is the inverse of the 
"tagged" cotangent CT K {x) a straightforward computation gives 



Recalling that 

i CT ^ - ■ 

and Clix) + kSI(x) 
mc 



x 2 + K 



hence 



dx 



X 2 + K 



-CT-\x). 



1 



S 2 (x) 



a 



1 the equations ([7]) and (fTT]) give 



f 



S 2 {cu + uq)' S 2 {cu + uq 

The given form of H is obtained by setting Wq = /o 



+ /o - ijikLq. 

niKLQ. 



Remark 4. The functions 7 and a = —rwy' in the non trivial case c 7^ take essentially three 
different forms that can be summarized by the three values k = 0, k = 1 and k = — 1. These 
three cases are described in the following table 







1 



7 



CU + Uq 



1 



tan(cu + no) 



tanh(cn + no) 



mc 



a 



(CU + Uq) 



mc 



a 



sin (cn + no) 



mc 



sinh (cn + no) 



Integrals Generated by Powers of an Operator 



7 



Lemma 3. For n > 2 the linear PDE © admits a (n+1)- dimensional linear space of solutions, 
if and only if the metric g has constant curvature K = mc. In particular the trivial case c = 
is possible only on a locally flat manifold. 

Proof. In components the equation © is 

ViVjG + mcgijG = d^G - T\fi k G + mcg^G = 0. 

Dividing by G and setting z k = dkG/G, since djkG/G = d k Zj + z k Zj, we get the first order 
PD-system in normal form (Pfaffian system) with n unknown depending on n variables 

di Zj — Zi Zj -\- r j j z k mcgij , 

whose integrability conditions are given by 

di(-zjzi + TflZk - mcgji) - dj(-ZiZ[ + T\z k - mega) = 0, (15) 

where the derivatives are computed by considering za as functions of (q 1 ). These conditions 
guarantee that the Pfaffian system is completely integrable, so that there exists a local solution 
for any choice of the values of the (zj) at a point qQ, i.e. a n-dimensional family of solutions 
parametrized by (a±, . . . , a n ) such that 

Moreover, since diZj = djZi any solution z\, . . . ,z n admits a potential Z such that Zi = d{Z = 
Sjln|G|. Thus, the solutions G of ([9]) depend also on an additional (n + l)-th multiplicative 
parameter. We expand (I15p and by introducing the Riemann tensor 

r>k o -pk o -pk , -ph -pk -phpk 

n lij — °i L jl ~ °j L ii "r 1 jiL ih — L il L jh 

and recalling that 

djgu - digji - Tjtfik + T a g jk = 0, 
we get 

R iij z k = mc(gjiZi - g a zj), 
that is, since they have to be satisfied for every z k , 

Rfij = mc( gj i5* - gu6*), 
which, by lowering the index k, are the constant curvature conditions 

Rhiij = K(gjig hi - gughj), K G R, 

and, therefore, K = mc. ■ 

Remark 5. Since g lj ' Hij = A, where A is the Laplace-Beltrami operator, the function G 
satisfies ([5]) only if it is an eigenfunction of the Laplace-Beltrami operator of (Q,g) with eigen- 
value —nmc, that is a solution of the Helmholtz equation with fixed energy. The condition is 
clearly not sufficient, because ([5]) must hold componentwise. 

Proposition 2. Let Li, 1 < i < 2n, k — 1 functionally independent first integrals of L in T*Q, 
let G, U m (G), V, H satisfying Theorem^] Then, for any m 6 N* 
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i) the k + 2 functions (H,U m (G), L, Li) are functionally independent if and only if 
{L,U m (G)}^0, 

ii) if L is a regular natural Hamiltonian on T*Q and G(q l ) is not constant, then {L,U m (G)} 

Proof, i) The rank of the Jacobian matrix of the (H,U rn (G), L, Li) w.r.t. the coordinates 
(u,p u , q l ,Pi) is equal to the rank of the square (k + 2) x (k + 2) matrix 



J 



1 a'L 


Pu 


dU m (G) 


dU m {G) 


du 


dp u 














V 






a dq a 


<) 


dU m (G) 


dU m (G) 


dq a 


dp b 


dL 


dL 


dq a 


dpb 


dL! 


8Li 


dq a 


dPb 



9L k 



8L k 



dq a 



dp b ) 



with \a\ + \b\ = k and where the indices a and b are chosen so that the rank of the k x k submatrix 
in the bottom-right corner, which we denote by is k. The determinant of J is given by 



det(J) = ( a'L 



dU m (G) dU m (G)\ , , » 

Because det(J^) 7^ by assumption, det(J) = if and only if the term between brackets is 
zero. This term is nothing but {H,U m (G)}N . Therefore, because {H,U m (G)} = 0, we have 
{H, U m (G)} Q = 0, but this is exactly equivalent to {L, U m (G)} Q = {L, U m (G)} = 0. 
ii) In the expression {L, U m (G)} the highest-degree term in p u is 

n 

-Y,d Pi Ld qi Gp™. 



If L is a regular (d Pi L ^ for all i) natural Hamiltonian and G is not constant then the term 
never vanishes, therefore {L, U m (G)} 7^ 0. ■ 

Corollary 1. If L is a regular natural Hamiltonian on T*Q and G is not constant, then the 
functions (H,U m (G), L, Li) of the previous proposition are all functionally independent. 

Corollary 2. Let L be a Liouville integrable Hamiltonian on T*Q with (Q,g) a (pseudo-)Rie- 
mannian manifold. If H, L, G satisfy Theorem^ then H is Liouville integrable with additional 
first integral U m (G) and, therefore, superintegrable. All additional first integrals of L are also 
additional first integrals of H . 

In the light of the previous results, in particular Theorem [21 the procedure for constructing 
extended Hamiltonians with at least one extra first integral U m (G) can be outlined as follows, 

1. Consider a constant-curvature (pseudo-)Riemannian manifold Q. 

2. Solve equation ([5]) for the functions G(q l ; a±, . . . , a n +i)- 

3. Solve equation ([6]) for the potential V and build L. 

4. Determine the extension through Proposition [2] and Remark [J] and fix the constant c and 
the integer m (with the constraint that their product is the curvature). 

5. Compute U m (G) to obtain the additional first integrals. 

In the following we provide some non trivial examples of the procedure outlined above. 
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3 Applications and examples 
Example 1: n = 1 

When Q is a one-dimensional Riemannian manifold there are no integrability conditions for the 
Hessian operator, moreover, after a rescaling of the coordinate, L can be always written as 

L = \pl + V{v). 

Hence, the geodesic part of the extended Hamiltonian H 

1 / 2 mc 

2 ^ + ^(cu + uq) 5 

corresponds to a Liouville metric on the manifold N x Q, with one Killing vector proportional 
to d v and constant curvature K = c 2 k. The conditions © and ([6]) of Theorem [1] become 

T-oGfw) = mcG(v), V'(v)G'(v) = -2meV(v)G(v). 
av A 

The solution of these differential equations are straightforwardly obtained by recalling that 

^-S K [x) = C K (x), ^-C K (x) = -kS k {x) 
ax ax 

and they are 

1 

G(v) oc S mc (v + u ), V(v) oc 



where oc means "is a constant multiple of". For c = 1 and k = 0, 1, —1, the manifold M = N xQ 
is respectively, the Euclidean plane, the sphere §2 and the pseudo-sphere Hfe. For c = — 1 and 
k = 0, 1, —1, the Minkowski plane, the deSitter and anti-deSitter manifolds, respectively. 

Remark 6. In [5] we obtained for the case n = 1 a dependence on m for the functions a, G 
and V, while here we have dependence on m for a only. However, a simple rescaling, that for 
this case is given by v — > — , shows that the two results are equivalent. This suggests that 
the most natural way to write the Tramblay-Turbiner-Winternitz like systems considered in 
HI HI CE21 [151 CE] is probably to shift the dependence on m from the potential function V 
of L to the metric factor a. A similar rescaling occurs for higher dimensional L. If L includes 
trigonometric functions of the (q % ), the rescaling makes evident the dihedral, polyhedral or, in 
general, discrete symmetries of L. 

Example 2: n = 2 

The flat case. If (Q,g) is a flat manifold with Cartesian coordinates (g 1 ,^ 2 ), a nonvanishing 
solution for G can be obtained by assuming c = (see Lemma [3]) and therefore 7 = — A(u + no) 
(LemmaH}. In this case, equation ([5]) becomes <9i<92G = and its solution is G = kQ+k\q l -\-k2q 2 , 
with ki not simultaneously zero. The compatible potential V obtained from equation ([6]) is 



V = niL, 



1 , k o \ 2 . ( 2 , k o 
? + 2k-J + { q + 2k 2 



+ F(k ig 2 -k 2 q 1 ), if hk 2 ^0, 



/ k \ 2 

V = mL U + -£) +F(q 2 ), if k 2 = 0. 
The two forms of V are equivalent up to rotations. 
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The pseudo-sphere. If (Q, g) = H2, the 2-dimensional pseudosphere of curvature K = — 1, 
with g 11 = 1, g 22 = 4(e ?? + e _,? ) _2 = cosh _2 (r/) in orthogonal coordinates (77, £), we have mc = — 1 
and from ([5]) 

G = (ai + a2e^ + a,3e~^)e~ v + (ci2e^ + a-^e~^ — a\)e n . 
For a2 = 03 = the integration of © gives for the potential 

no 



V 



cosh 2 (77) 



In this case V is in Stackel form V = g n fi{q l ) and compatible with the separation of variables 
of L, which is therefore integrable with one quadratic first integral H\ at least. The extended 
Hamiltonian H is therefore always superintegrable with the four first integrals H, L, Hi and 
U m (G) for mc = —1. We write explicitly U m (G), m = 1,2, relatively to the extension 

1 2 m 2 m 

H = T,Pu 2~ L ' ^ = ' 

as follows 

U{G) = -2smh(r,) Pu + 2^^-pi, 

u 

U\G) = -2sinh(^ + 8 C -^l PuPl - 16^Ml. 

u u z 

The sphere. If (Q,g) = § 2 , the 2-dimensional sphere of curvature K = 1, with g 11 = 1, 
g 22 = sin -2 9) in standard spherical coordinates (9, (p), we have mc = 1 and from © 

G = {a\ sin <fi + 0,2 cos (f) sin 9 + 03 cos 9. 

For 03 = 0, the integration of © gives for the potential 

V = 7t-r F ((01 c °s (p — (12 sin^) tan 9) , 

cos z 9 

which in general is not in Stackel form and therefore we have an example of a non separable L 
which admits an extension (again not separable) with a polynomial additional first integral. For 
m = 1 and F = (a\ cos <j) — (12 sin <fi) tan 9, for example, we have with V of above 

U (G) = (ai sin (j) + d2 cos <p) \ Pu sin# H — P $ cos 9) H ^P<a(oi cos (f> — 02 sin<^>). 

\ u I it sin 



Example 3: n = 3 

Let us consider as (Q,g) the sphere § 3 with coordinates (q 1 = m q 2 = £1, g 3 = £2) where 
< 77 < 7r/2 and < £j < 27T and the parameterization in ]R 4 being given by 

x = cos£i sin#, y = sin£i sin#, z = cos£2 cos#, i = cos£2Cos#. 

These coordinates are known as Hopf coordinates, the non null components of the covariant 
metric tensor are gn = 1, 522 = sin 2 77 and 533 = cos 2 7/ so that the curvature is -K" = 1. The 
surfaces r/ = const are flat tori spanned by the rotations which are Killing vectors of the 
manifold. These coordinates correspond to a cylindrical rotational separable system, they are 
associated with a Killing 2-tensor K which is described in appendix in [B] and which provides 
a quadratic first integral of the geodesies Hi = ^K^ P i P j. Therefore, the geodesic Hamiltonian 
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G = \g ll Pi of S3 admits the following four independent quadratic in the momenta first integrals 
G, Hi, p^, p\ and is Liouville integrable. By applying our procedure we obtain from ([SD 

G = (03 sin£i + 04 cos£i) sin 77 + (a± sin £2 + 0-2 cos £2) cos rj, 

with aj constants. After setting 02 = 03 = 04 = one obtains easily from ([6]) 

sin 77 V cos £2 

Remarkably, this potential can be compatible with separation of variables, then making G+V an 
integrable system, if it is in Stackel form V = g n fi(q l ) as happens for example with F = F(£i). 
However, in general it is not, as for 



sin£ 



cos £ 



-tan?/. (16) 



In both cases our procedure provides independent first integrals U m (G) for the extended Hamil- 
tonians H described in the previous sections, with the prescription mc =1. In the case of 
F = F(£i), for example, H is a superintegrable Hamiltonian with, at least, four other inde- 
pendent first integrals. In the case (|16fl . H admits the two other independent first integrals L 
and U m (G) at least. For V with F given bv (|16p . we have for example 

TTtn\ ■ c , 1 ( cos & ■ • t 

U (G) = p u cos r/ sin £ 2 + - P3 Pi sin rj sin £ 2 

U \ COST] 



4 Conclusions 

We have proved that the method developed in [5] for building polynomial additional first inte- 
grals of arbitrarily high order of extended two-dimensional Hamiltonians can be generalized to 
similar extensions of n-dimensional Riemannian or pseudo-Riemannian manifolds. As a conse- 
quence of natural integr ability conditions, it is necessary that the curvature of these manifolds 
is constant. In examples for n = 1,2,3 we improve the results obtained earlier for the one- 
dimensional manifolds and show how the procedure works in the two and three dimensional 
ones. Future directions of research will be towards the characterisation of more general proce- 
dures of extension of Hamiltonian systems and the search for more general expressions of the 
operator U and of the function G. 
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